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We study the second-order Raman process of mono- and few-layer MoTe 2 , by combining ab initio 
density functional perturbation calculations with experimental Raman spectroscopy using 532, 633 
and 785 nm excitation lasers. The calculated electronic band structure and the density of states 
show that the resonance Raman process occurs at the M point in the Brillouin zone, where a strong 
optical absorption occurs due to a logarithmic Van-Hove singularity of electronic density of states. 

Double resonance Raman process with inter-valley electron-phonon coupling connects two of the 
three inequivalent M points in the Brillouin zone, giving rise to second-order Raman peaks due to 
the M point phonons. The calculated vibrational frequencies of the second-order Raman spectra 
agree with the observed laser-energy dependent Raman shifts in the experiment. 

PACS numbers: 33.40.+f, 78.30.Hv, 63.20.kd, 63.22.Np, 63.22.-m, 31.15.A- 


I. INTRODUCTION 


Transition metal dichalcogenides (TMDs) have been 
intensively studied over past decades, and recently 
attracted tremendous attention for a number of 
applications^ - — , since TMDs have a finite energy band 
gap, which is more advantageous than graphene. Among 
the TMDs, MoTe2 has a stronger spin-orbit coupling 
than other Mo related TMD materials (M0X2, X = S, 
Se), which leads to a longer decoherence time both for 
valley and spin freedom^ and thus MoTe2 is considered 
to be used for ideal valleytronic devices^. Unlike other 
TMDs, MoTe2 has unique structural properties, such 
as a thermally-induced structural phase transition at 
around 900°C from the <a-phase diamagnetic semiconduc¬ 
tor to the /3-phase paramagnetic metal&£. MoTe2 has the 
smallest direct energy band gap 1.1 eV) at the K point 
(hexagonal corner) in the Brillouin zone (BZ) among all 
semiconducting TMDs, the band gap is similar to the 
indirect gap (~ 1.1 eV) of silico n 10 ! 11 . MoTe2 is thus 
suitable for nanoelectronics devices^ and for studying 
exciton effects at the K points^. Raman spectroscopy 
is a useful tool for characterizing TMDs since we get 
many information^ - — such as number of layers, elec¬ 
tronic and phonon properties of TMDs. Recent Raman 
spectroscopy studies of monolayer and few-layer MoTe2 
using visible light (1.9 ~ 2.5 eV) show that there ex¬ 
ists some unassigned Raman spectra which can not be 
assigned to a first-order Raman process**^, and the as¬ 
signment of these Raman spectra is the subject of the 


present paper. 

Raman spectra of TMDs have been reported for sev¬ 
eral decadesii - — . The first-order Raman process in 
TMDs has been widely investigated and is now well 
understood^ ^ i ^^ i ^ r— . However, the study of second- 
order Raman spectra did not give consistent assignments 
for different TMD a 18 i 25 i 26 i 35 — . For example, Stacy et 
al 26 and many other group s 15 ! 18 ! 29 ! 38- — assigned some 
of the Raman peaks of 2H-M0S2 to either combination 
or overtone second-order phonon modes occurring at the 
M point (center of hexagonal edge) in the Brillouin zone. 
Terrones et al£l also assigned phonon modes in few- 
layer WSe2 according to the M point phonon. On the 
other hand, Sourisseau et al 27 ' 28 ascribed the second- 
order Raman process of 2H-WS2 to the K point phonons. 
Berkdemir et considered that electrons are scattered 
by the M point phonons between the K point and I point 
(one point along TK direction), and that this process 
is responsible for the double resonant Raman process in 
few-layer WS2. Furthermore, Sekine et alS^ found a dis¬ 
persive Raman peak in 2H-M0S2 and interpreted this 
peak in terms of a two-phonon Raman process at the 
K point. The different assignments for both few-layer 
or bulk M0S2 and WS2, either at the M or K points, 
have not been well investigated, since most of these as¬ 
signments to second-order Raman peaks have been made 
merely by comparing phonon frequencies with reference 
to some inelastic neutron scattering result s 19 ! 28 . Thus 
a detailed analysis with use of double resonance Ra¬ 
man theory is needed, especially for understanding both 
mono- and few-layer MoTe2 that are discussed in this 
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paper. 

In the present paper, we mainly study the second-order 
Raman spectra of monolayer MoTe 2 , based on ab initio 
electronic band calculations, density functional pertur¬ 
bation theory and associated experimental Raman spec¬ 
troscopy. From our electronic band structure calcula¬ 
tion, we have found a saddle point of the energy disper¬ 
sion at the M points in the BZ, which gives rise to two- 
dimensional Van Hove singularity (VHS) in the electronic 
density of states. The energy separation between two 
VHSs (2.30 and 2.07 eV) are matched to the 2.33 eV (532 
nm) and 1.96 eV (633 nm) laser excitations in our ex¬ 
periment in which resonant optical absorption occurs at 
the M points. Consequently, inter-valley electron-phonon 
scattering process occurs between two of the three in¬ 
equivalent M (M, M', M") points, giving rise to a double 
resonance Raman process. We propose that this double 
resonance process occurring at the M points is essential 
for the interpretation of the observed second-order Ra¬ 
man peaks of MoTe 2 - The experimental second-order 
Raman frequencies as a function of laser excitation en¬ 
ergy confirm our interpretation of the double resonance 
Raman spectra. 

Organization of the paper is as follows. In II and III, 
we briefly show computational details and experimental 
methods, respectively. In IV, both experimental and the¬ 
oretical results for mono- and few-layer MoTe 2 Raman 
spectra are given. We assign second-order Raman peaks 
by comparing the calculated results with the observed 
Raman spectra. In V summary of the paper is given. 
Further information such as layer-number dependence is 
discussed in Supplemental Material. 


II. COMPUTATIONAL DETAILS 

We performed the electronic and phonon energy dis¬ 
persion calculations on monolayer and few-layer MoTe 2 
by using first-principles density functional theory within 
the local density approximation (LDA) as implemented 
in the Quantum-Espresso code^. The structures of 
monolayer and few-layer MoTe 2 that are used in the 
present experiment are a-MoTe^. The A A' and A A'A 
stacking geometries are used in the computational study 
for bi-layer and tri-layer MoTe 2 , respectively, since these 
stacking geometries are more stable than the other 
possibilities^. AA' stacking refers to the geometry which 
has a Mo (Te) atom of one layer on top of a Te (Mo) atom 
of the other layer. The few-layer MoTe 2 are separated by 
25 A from one another in a unit cell of the calculation 
to eliminate the inter-few-layer interaction. We used pro¬ 
jector augmented-wave (PAW) pseudopotentials^ with a 
plane-wave cutoff energy of 65 Ry to describe the inter¬ 
action between electrons and ions. The spin-orbit split 
electronic band structures were calculated by relativistic 
pseudopotentials derived from an atomic Dirac-like equa¬ 
tion. The atomic coordinates are relaxed until the atomic 
force is less than 10 -5 Ry/Bohr. The lattice parame¬ 


ters used in this paper are listed in Table SI and agree 
well with the experiment and other calculation s 32 ! 45 ! 46 . 
The Monkhorst-Pack scheme^ is used to sample the Bril- 
louin zone (BZ) over a 17x17x1 £;-mesh for MoTe 2 - The 
phonon energy dispersion relations of MoTe 2 are calcu¬ 
lated based on density functional perturbation theory^. 
In order to visualize the vibrational modes at the M 
point, we use a 2x2x1 supercell in which the M point 
is zone-folded to the T point in the folded BZ. 

The optical absorption spectrum is calculated by the 
real (e') and imaginary (e") parts of the dielectric func¬ 
tion as a function of photon energy, respectively, based on 
the PAW methodology^ and the conventional Kramers- 
Kronig transformation. The absorption coefficient a is 
described by a = 47ncE 7 /(ftc), where E 7 is the incident 
laser excitation energy, h the Planck constant, c the speed 
of light, and k is the extinction coefficient^, that is, k = 

\j (V e' 2 + e" 2 — e')/2. 


III. EXPERIMENTAL METHOD 

Bulk crystals of <a-MoTe 2 were prepared through a 
chemical vapor transport method^. Atomically thin 
crystals of MoTe 2 were mechanically exfoliated from the 
bulk crystals onto silicon substrates with a 90 nm-thick 
oxide layer on top. The thicknesses of the atomic crystals 
were identified using optical microscopy and first-order 
Raman spectroscopy^. Raman spectroscopy measure¬ 
ments for monolayer MoTe 2 were performed using 532, 
633 and 785 nm excitation lasers for discussing the ob¬ 
served phonon dispersion due to double resonance Ra¬ 
man theory 7 . Raman spectroscopy studies of few-layer 
MoTe 2 under 532 nm laser excitation were also performed 
and the results are discussed in the supporting materials. 
The grating sizes were 1200, 1800 and 2400 lines/mm 
for the 785, 633 and 532 nm laser excitation measure¬ 
ments, respectively. The magnification of the objective 
lens was lOOx. The accumulation times were 150-200 
seconds. The laser power was kept below 0.1 mW, 0.2 
mW and 3.0 mW for the 532, 633 and 785 nm excita¬ 
tion lasers, respectively, in order to avoid damage to the 
samples. All measurements were performed at room tem¬ 
perature in the backscattering configuration. Since this 
paper gives greater emphasis to a theoretical understand¬ 
ing of the 2nd-order Raman scattering processes, we do 
not present many experimental details. However, more 
details about the characterization of the number of lay¬ 
ers in our MoTe 2 samples can be obtained in one of the 
previous works on MoTe 2 synthesis and characterization, 
based on the first-order Raman spectral. 


IV. RESULTS AND DISCUSSION 

Figure [lja) shows Raman spectra of mono- to tri¬ 
layer MoTe 2 under the 532 nm (2.33 eV) laser excitation. 
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(a) (b) 



Raman shift (cm 1 ) Raman shift (cm 1 ) 

Figure 1: (Color online) (a) Raman spectra of mono-, bi-, and tri-layer MoTe 2 under 2.33 eV (532 nm) laser excitation energy, 
showing the small intensity peaks uji (i = 1, 2, ..., 6) and an additional calibration peak from Si besides the typical Raman 
peaks of Ai^, E\ g and E\ g that are frequently observed in MoTe^^. (b) Raman spectra of monolayer MoTe 2 under 2.33 eV 
(532 nm), 1.96 eV (633 nm), and 1.58 eV (785 nm) laser excitation energies. 


The Raman spectra show three strong peaks, including 
the in-plane E\ g mode and the out-of-plane Ai^ mode 
for monolayer and few-layer MoTe 2 , as well as the E\ g 
mode and the E\ g mode for bi- and tri-layer MoTe 2 , that 
were observed previously^ ! 16 ! 52 . These spectra are as¬ 
signed to the first-order Raman spectra of the T point 
phononsi^. However, we observe additional peaks with 
relatively small intensities, which can not be assigned as 
first-order Raman spectral. In Fig. Oja), we denote the 
weak Raman spectra as u)i (i = 1, 2, ..., 6). 

Since the second-order Raman spectra are dispersive 
as a function of excitation energy 7 , we performed Ra¬ 
man spectroscopy measurements using different excita¬ 
tion laser energies. In Fig. [ljb) we show Raman spectra 
of monolayer MoTe 2 under 2.33 eV (532 nm), 1.96 eV 
(633 nm) and 1.58 eV (785 nm) laser excitation energies. 
Here the Raman intensities are normalized using the non¬ 
resonant Si peak intensity at 520 cm -1 (not shown in 
Fig.[ljb)). We find that the peak positions of each uji ei¬ 
ther upshift or downshift by up to several cm -1 for differ¬ 
ent laser excitation energies, suggesting that these peaks 
are associated with a second-order Raman process 7 . 

Considering that the second-order process is usually 
connected to the electronic properties by double reso¬ 
nance Raman theory 7 , we calculate the electronic band 
structure, the density of states (DOS) and the optical ab¬ 
sorption spectra of monolayer MoTe 2 as shown in Fig.[2j 
Our calculated energy dispersion reproduces the main ob¬ 
served feature s 11 ! 32 ! 54 such as the direct electronic band 
gap (1.00 eV) and the spin-orbit splitting of the valence 
band (^230 meV) at the K point which are close to 1.10 
eV and 250 meV, respectively, obtained from the op¬ 
tical absorption measurements^. It is noted that the 
energy band gap is underestimated by the local den¬ 
sity functional method. In order to fit the experimental 
data, the conduction bands are upshifted by 0.10 eV, 
so are the DOS and optical absorption spectrum. In 


Fig. [2ja), we show the calculated electronic energy band 
and DOS. From the parabolic bands around the K point, 
two-dimensional (2D) DOS gives constant values, while 
at the M point we have found saddle points of the energy 
dispersion in the electronic band, the corresponding DOS 
gives logarithmic log|E - E$\ divergence as is known as 
Van Hove Singularity (VHS) in 2D systems. 

In Fig.[2jb), we show the calculated real (e') and imag¬ 
inary (e") parts of the dielectric constant and the optical 
absorption coefficient aasa function of photon energy. 
e" has a nonzero value from the photon energy E 7 = 1.10 
eV which corresponds to the energy gap, and a peak at 
E 7 ~ 1.80 eV. A strong optical absorption at 1.80 eV 
from our calculation agrees well with the large peak at 
around 1.83 eV from the optical absorption measurement 
by Ruppert et al*ii on monolayer MoTe 2 - Moreover, the 
electronic excitation energy of 2.30 eV (2.07 eV) from the 
third (top) valence band to the lowest (third lowest) con¬ 
duction band at the M point matches very well with the 
laser energy 2.33 eV (1.96 eV) used in our experiment, as 
indicated by pink (orange) arrows in the band structure 
of Fig. mja). Thus resonant optical absorption at these 
laser lines is expected. 

In Fig. [2jc), we plot equi-energy contour for the con¬ 
duction band energies E^ an ^ which correspond to the 
laser excitation energies E 7 . E b an d = 1.52 and 1.57 eV for 
cl and 1.72 eV for c2 correspond to the E 7 = 1.58, 2.33 
and 1.96 eV, respectively, as indicated by blue, pink, and 
orange arrows in Fig. [2fa). The conduction bands cl, c2 
and the valence bands which have the same spir&£2 are 
highlighted in red solid lines in the band structure. The 
equi-energy contour lines Eg 1 ^ and Eg 2 nd in Fig. Elc) 
both show the saddle point shape of the energy dispersion 
in the vicinity of the M point. That is, the energy band 
E ^nd (E band ) starting from the M point monotonically 
decreases (increases) along the M-K direction, whereas 
this band energy increases (decreases) along the M-T di- 
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Figure 2: (Color online) (a) Electronic energy band structure E(k) and the density of states (DOS), (b) the real (e x ), and 
imaginary (e 7/ ) parts of the dielectric function and the optical absorption coefficient a as a function of photon energy E 7 , and 
(c) an equi-energy contour plot of the two conduction bands (cl and c2) with the same spin for monolayer MoTe 2 . Energy 
bands are split by spin-orbit interaction, the conduction bands cl, c2 and the valence bands which have the same spin^^ 
are highlighted in red solid lines in (a). When we consider energy separation between two Van-Hove singularities (VHSs) 
of the DOS, the peak at 1.80 eV in the absorption spectrum corresponds to optical transition at the M point (black double 
headed arrow in (a)), (c) The equi-energy contour plot for VHS E^ nd = 1.52 and 1.57 eV (left) and E \? and = 1.72 eV (right). 
All equi-energy contours appear around the M point. The E band energies of 1.52, 1.57 and 1.72 eV (horizontal lines in (a)) 
correspond to the experimental laser energies (1.58, 2.33 and 1.96 eV) in blue, pink and orange colors. 


rection. When we decrease laser excitation energies from 
2.33 to 1.58 eV, the contour lines in Fig.[2fc) left change 
the shapes around the M point, suggesting that the area 
inside of the contour line for each M point decreases with 
increasing laser excitation energy. This situation is op¬ 
posite to the case of graphene in which the equi-energy 
contour is almost circular and the area increases with 
increasing energy. 

If the laser excitation energy matches the optical tran¬ 
sition energy at the M point, we expect that many photo- 
excited electrons at the M point are scattered to the in¬ 
equivalent M' or M" points by an inter-valley resonant 
electron-phonon interaction, as shown in Fig. [3^a). The 
inter-valley electron-phonon process, together with the 
electron-photon process, gives rise to a double resonance 


Raman process. Selecting one process from three possible 
inter-valley scattering of photo-excited carriers (Mi — >M' 
or Mi —or M'i —^M"), we show the schematic view 
of the double resonance process in Fig.[3jb). In Fig.[3jb), 
a photo-excited electron will be scattered from M to M' 
by emitting a phonon with energy huj p i and momentum 
—g, which conserves the energy and momentum of both 
the electron and the phonon. Then the electron is scat¬ 
tered back to M by emitting the second phonon with 
energy huo p 2 and momentum +q. As analyzed in the 
Fig. Ha), the momentum q of the phonon can only have a 

magnitude of |MM'| or |MM^| or | M / M / ^ |. When we mea¬ 
sure the value of |q| from the T point, the q corresponds 

to rM ;>/ or TKM' or TM, respectively. Since rKM 7 = 
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Figure 3: (Color online) (a) In the double-resonance pro¬ 
cess, photon-excited electrons are scattered inelastically by 
emitting phonons and make inter-valley (i.e., M< —>-M 7 , 
M' < —>-M 77 , and M< —>M 77 ) transitions, as shown by the solid 
double-end arrows in the first BZ. The related phonon vectors 
q measured from the T point are also given by the six dashed 
arrows, (b) The schematics of the double resonance process 
(electron-photon and electron-phonon processes). This pro¬ 
cess generates two phonons (huj p i, -q) and {frw p 2 , q)- 


TM' + F? = TM', all phonon modes related to double 
resonance process correspond to phonon modes at the M 
point. 

In order to specify the phonons contributing to the 
second-order Raman process, we analyzed possible inter¬ 
valley phonon vectors qMM' f° r the three laser energies 
and calculated the phonon spectra of monolayer MoTe 2 , 
as shown in Fig. 0J For the laser with an energy of 2.33 
eV, we analyze possible qMM' wave vectors between two 
elliptical equi-energy contour lines (red lines) in the left 
panel of Fig. SJa). The dark area for g MM / shows that 
the most probable phonon vectors appear at the M point. 
Other phonon vectors, such as qM'M" and qMM", are sim- 
ilar to gMM' (not shown). For the 1.96 and 1.58 eV exci¬ 
tation lasers, as already seen in Fig.[2fc), the equi-energy 
contours have a larger area than that for 2.33 eV. Then 
we expect that qMM' f° r E 7 = 1.96 and 1.58 eV arises 
from the vicinity of the Brillouin zone near the M point, 
as shown in the central and right panels of Fig. |U^a). 
The relatively dark region indicates that: (1) the qMM' 
for E 7 = 1.96 eV appear both at the M point and in the 
vicinity of the M point (approximately with lengths of 
|MK) along the MK direction and (2) the qMM' for E 7 
= 1.58 eV appear in the vicinity of the M point along the 
M 7 M 77 (or MM 77 ) direction approximately with lengths of 
10%M / M // (or 6.5%MM // ) away from the M 77 point. 

In Fig.lUb), we show the calculated phonon dispersion 
relations of monolayer MoTe 2 - The symmetry for each 
optical phonon band?^, such as E \ g and Ei^, is labeled in 
Fig.lHb). Here LA, TA and ZA refer to the longitudinal, 


transverse and out-of-plane acoustic modes, respectively. 
We assign the second-order Raman modes at the M point 
for the laser energy of 2.33 eV. The strategy for the as¬ 
signment is to choose any two phonon modes at the M 
point as one overtone or (difference) combination mode 
and to list in Table |T] those that have frequency values 
within the experimental 07 =L 6 cm -1 . These assignments 
are marked in Fig. S^b) by using color arrows. The red, 
green, blue, purple, magenta, and cyan arrows represent 
overtone or combination Raman modes of 07 (i = 1, 2, ..., 
6 ), respectively. Downward arrows with an a '” super¬ 
script represent a two-phonon process with one phonon 
annihilated and one phonon created to generate ” differ¬ 
ence combination mode”, while upward arrows with the 
” +” or ” 2 ” superscripts are used for combination or over¬ 
tone modes in which two phonons are created. For exam¬ 
ple, uj 1 is assigned to either the cjf (= E^(M) - LA(M)) 
or co\ (= 2TA(M)), cj 2 to uFT (= F^(M) - TA(M)), u ; 3 
to either cj§ (= 2 LA(M)) or (= E lp (M) + TA(M)), 
0 J 4 to c J 4 " (= Ai^(M) +LA(M)), UJ 5 to either the (= 
2 Ai^(M)) or CJ 5 - (= E^(M) + TA(M)) and ujq to the uJq 
(= E^(M) + LA(M)) mode. Using the same phonon as¬ 
signment, we obtain the second-order Raman frequencies 
for E 7 = 1.96 eV both at the M point and near the M 
point with a length of |MK and for 1.58 eV at one point 

away from the M point with 10%M / M // . 

In Table [U we list phonon frequencies of the first-order 
(upper) and the second-order (lower) Raman spectra for 
both the experiment and the calculation and also the 
phonon assignments. The first-order experimental Ra¬ 
man spectra (Ai p and E^) do not change the frequency 
for the laser energies (2.33 and 1.96 eV) within the exper¬ 
imental error (1 cm -1 ) and agree with the calculations. 
The small frequency downshift of the first order spectra 
(1.7 cm -1 at most) from 2.33 eV to 1.58 eV laser excita¬ 
tion is checked to be not from the laser-induced thermal 
effect as discussed in the Supplemental Material. Al¬ 
though we did not specify the origin of the downshift 
precisely, the origin of the change for the first order can 
be neglected for simplicity. The experimental second- 
order peak positions of 07 are either upshifted or down¬ 
shifted with the laser energy E 7 by a small magnitude 
of about 1-3 cm -1 . It is pointed out that the frequency 
of 0 J 1 is upshifted from 2.33 eV to 1.58 eV even though 
the both the first-order spectra are downshifted. This 
can be understood by double resonance Raman theory 
in which the resonance phonon frequency changes along 
the phonon dispersion near the M point. In the right of 
Table [H the second-order phonon frequencies for 2.33 and 
1.58 eV laser excitations are obtained at qMM' = TM (at 
M) and ^10%MM / (denoted by ”off M”), respectively. 
For E 7 = 1.96 eV, the phonon assignments are obtained 

either at qMM' = Tivl (at M) or ~ (off M). 

These assignments can be further confirmed by the ob¬ 
served relative Raman intensity and the phonon disper¬ 
sion relation. First, in Fig.[]Jb), a larger relative intensity 
of Raman spectra 07 can be found at E 7 = 2.33 eV than 
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2.33(1.96) eV 1.58 eV 1.96 eV 



Figure 4: (Color online) (a) All possible phonon vectors q M M' coming from the double resonance Raman process due to the 
laser energies of 2.33, 1.96 and 1.58 eV. The <f M M' is equal to either FM (for the laser energy E 7 = 2.33 and 1.96 eV) or 
~ |MK (for E 7 = 1.96 eV) or ~10%MM / (for E 7 — 1.58 eV), as indicated by an arrow away from the F point, (b) The phonon 
dispersion relations of a monolayer MoTe 2 - The second-order phonon modes are analyzed at the M 7 point for the laser energy 
E 7 = 2.33 eV, and then at two k points marked by two vertical dashed lines which are shifted leftward and rightward from 
the M point for the laser energy E 7 — 1.58 eV and 1.96 eV, respectively. The red, green, blue, purple, magenta, and cyan 
vertical arrows are used for modes tc?i, CJ 2 , CJ 3 , CJ 4 , CJ 5 , and u> 6 , respectively. The superscripts ”+”, and ”2” for uj refer to 
possible ”combination modes” (i.e. = Ai g (M) + LA(M), ujq — E^(M) + LA(M)), ’’difference modes” (i.e. cjj - = E^(M) 

— LA(M), UJ 2 = E^(M) — TA(M)) and ’’overtones” (i.e. uj\ — 2TA(M), — 2LA(M)), respectively. 


at 1.58 eV. This supports the assignments at the M point 
and at q MM ' r ^10%MM / for E 7 = 2.33 eV and 1.58 eV, 
respectively, since phonon density of states (PhDOS) at 
the M point is larger than that at ^10%MM / as 

shown in Fig. HJ^b). Similarly, for E 7 = 1.96 eV, a rela¬ 
tively large (small) intensity of the Raman spectra 002 , ^3 
and ujq (c^q, CJ 4 and CJ 5 ) in Fig. [Tfb) may arise from the 


large (small) PhDOS at the M point (off the M point). It 
is therefore reasonable to assign the CJ 2 , ^3 and (cji, 
CJ 4 and CJ 5 ) at the M point (off the M point), as under¬ 
lined in Table [H Second, it is noted that the experimental 
spectra of CJ 4 (at the three laser energies) are too weak 
to determine the frequency values. The uj 4 is assigned 
as a combination mode related to the Ai p (M) as listed 
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Table I: Phonon frequencies of the first- and second-order Raman spectra of monolayer MoTe 2 are listed for three different 
laser energies E 7 . A comparison for first and second-order Raman spectra is given between the experimentally unassigned 
modes uji (i = 1, 2, 6 ) and the calculated phonon frequencies both at the M point and near the M point denoted by (off M) 

whose q lies around |MK and 10%MM' for E 7 = 1.96 eV and 1.58 eV, respectively. The consistent dispersion behavior for the 
second-order Raman spectra between the experiment and the theory are underlined. The frequency is given in units of cm -1 . 
NA: not available. 



Experiment 

Calculation 

Phonon assignments 

E 7 

2.33 eV 

1.96 eV 

1.58 eV 

2.33 eV 

1.96 eV 

1.58 eV 


El g 

NA 

NA 

NA 

116.9 



first-order 

Alg 

171.1 

171.3 

170.0 

172.2 



first-order 

K 

235.8 

236.2 

234.1 

236.9 



first-order 



(at M) 

(at M/off M) 

(off M) 


Ul 

138.2 

140.7 

139.0 

137.4 

137.4/146.6 

140.1 

2TA(M) 





132.8 

132.8/133.1 

133.2 

or E 2 9 (M)-LA(M) 

UJ2 

183.7 

183.6 

181.5 

178.3 

178.3/177.2 

177.8 

Ei,(M)-TA(M) 


205.5 

205.3 

205.9 

197.9 

197.9/208.2 

200.5 

Ei 9 (M)+TA(M) 





202.7 

202.7/196.0 

200.4 

or 2LA(M) 

UJ4 

-264.0 

260-270 

260-270 

263.0 

263.0/256.6 

261.0 

A l 9 (M) + LA(M) 

^5 

329.0 

333.2 

-328.0 

323.3 

323.3/317.2 

321.6 

2Ai 9 (M) 





315.7 

315.7/323.8 

317.9 

or E^ 9 (M)+TA(M) 


347.4 

348.4 

346.5 

348.3 

348.3/348.5 

348.0 

E^ 9 (M)+LA(M) 


in Table |T| Then it is interesting to discuss why the uo 4 
becomes weaker and broader when the first-order A ig( r ) 
mode becomes stronger by changing the laser excitation 
energy from 2.33 eV to 1.96 eV/1.58 eV. We expect that 
the photo-excited electron may choose either an intra¬ 
valley scattering process (Ai g (T)) or an inter-valley scat¬ 
tering process (Ai p (M)), so that Ai g (M) competes with 
Aip(T), resulting in the alternative Raman intensity of 
Aip(T) and CJ4 in the experiment. The CJ4 (Ai g (M)) has 
no VHS of PhDOS compared with all the other cj’s at 
the M point as seen from Fig. Hfb) which is relevant to 
the weak intensity of CJ 4 , too. 

In Fig. [5j we show the normal modes of the phonons 
at the M point in which the phonon vibration ampli¬ 
tude is visualized in the 2 x 2 supercell. As seen from 
Fig. |4jb), the two groups of optical phonon modes (E^ 
and (Ei p or Ai p )) are well separated by about 50 cm -1 . 
The visualization of these M point phonon modes shows 
that the LO and TO E 2 ^(M) modes (230^300 cm -1 ) and 
the LA(M) and TA(M) modes (60^100 cm -1 ) are bond¬ 
stretching modes of the Mo-Te bond. The LO Ai g (M) 
and TO Ei p (M) modes (120^180 cm -1 ) arises from the 
relative motion of two Te atoms occupying the same sub¬ 
lattice. The six vibrational modes in Fig. [5] are Raman 
active, since the volume or the area of the supercell un¬ 
der vibration gives non-zero deformation potential. Thus 
all the modes can contribute to the 2 nd-order overtone, 
combination and difference modes. Furthermore, the vi¬ 
brational magnitude of the TA modes is no less than that 


of the LA modes, suggesting that the TA modes may in¬ 
duce as large deformation potential as the LA modes in 
the electron-phonon interaction. Usually, the TA modes 
around the T point do not contribute to the electron- 
phonon interaction. Nevertheless, this seems not to be 
the case at the M point. 

So far we have only considered the electron-involved 
double resonance process for MoTe 2 and did not dis¬ 
cuss the hole-involved double resonance process. As 
Venezuela et alJ& pointed out, the enhancement by the 
double resonance is strong in graphene for one-electron 
and one-hole scattering, since the energy denominators 
of the Raman intensity formula becomes zero for the 
symmetric it and 7r* bands of graphene (triple reso¬ 
nance). In the case of TMDs, however, a similar effect 
is not expected since the conduction band and the va¬ 
lence band are not symmetric around the Fermi energy^. 
Since we don’t have a program to calculate electron- 
phonon matrix element by first-principles calculations, 
we can’t calculate Raman intensity for each double reso¬ 
nance process^ which will be a future work. 

Finally, we briefly comment on the exciton effect on the 
double resonance Raman spectra of MoTe 2 - Since the ex- 
citonic states appear at the direct energy gap at the K 
point of the B Z 56 i 57 , we expect that the excitonic energy 
is around 1.1 eV for MoTe 2 , which is much smaller than 
the three laser energies we used (1.58 - 2.33 eV). As far as 
we consider only MoTe 2 , we can safely neglect the exciton 
effect at the K point. However, in the case of other TMDs 
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Figure 5: (Color online) Visualization of the possible Raman- 
active modes at the M point for the second-order process sug¬ 
gested in Figure 4. Both the top and side views are given 
and red arrows represent the vibrational directions and am¬ 
plitudes of the motion of individual atoms. 


Table II: Laser energies that are used in the Raman spectra of 
TMDs. Ec,(K), E^(M) are, respectively, the optical transition 
energies at the K and M point, which are obtained by LDA 
calculations. Detail is discussed in the text. 


TMDs 

laser (eV) 

E 9 (K) E 9 (M) 

MoSsr^ 

1.83-2.41 

1.7 

>3.0 


1.83-2.41 

1.9 

>3.0 

WSe^^ 

2.41 

1.7 

>3.0 

MoTe 2 this work 

1.96, 2.33 

1.1 

-1.8 


(MoS^, WS ^W 1 ! 58 , WSe 2 36l3T ) energy gaps at the K 
point (E^(K)) lie in the range of 1.7-1.9 eV as shown in 
Table HI The laser energies used in the experiments lie 
in a range from 1.83 to 2.41 eV. Thus the exciton effect 
needs to be considered for some laser energies. Neverthe¬ 


less, when we discuss the M-point phonon of the other 
TMDs, we expect that the second-order Raman spectra 
do not appear, since the optical transition energies at the 
M point is larger than 3.0 eV. It would be interesting if 
ultra-violet laser could be used for detecting the M-point 
phonon modes for other TMDs, which will be a future 
work. In the case of MoTe 2 , on the other hand, since 
Efl(M) 1.8 eV, we expect the double resonance Raman 
process for conventional laser energies. Although we did 
not consider the exciton effect at the M point explicitly 
as K.F. Mak et alJ& did in graphene, we believe that the 
present discussion of the double resonance Raman effect 
should work well, because only optical transition matrix 
elements that are common in the resonant Raman pro¬ 
cess are enhanced by the exciton effect, and the relative 
Raman intensity thus obtained do not change much for 
the M point 


V. CONCLUSION 


In summary, we have investigated the second-order Ra¬ 
man spectra of MoTe 2 systems, based on an ab initio 
density functional calculation, and compared with the 
experimental Raman spectroscopy. The calculated elec¬ 
tronic bands, the Van Hove singularities in the density of 
states, and the absorption spectra indicate a strong opti¬ 
cal absorption at the M point in the BZ of MoTe 2 - Using 
double resonance Raman theory, combination/difference 
or overtone mode of phonons at the M point are assigned 
to identify the experimental Raman modes of uji (i = 1, 
2, ..., 6). This study facilitates a better understanding 
of the electron-phonon interaction and the second-order 
Raman process in TMD systems. 
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